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SUMMARY 


The major effort of the ALPEX project at the University of Wisconsin 
Space Science and Engineering Center was software development of TOVS process- 
ing on the Man-computer Interactive Data Access System (McIDAS) initially 
based on the IBM 4341 and now on the IBM 4381. This work included: (1) the 

addition of a high resolution topography, (2) formulation and implementation 
of a new surface skin temperature calculation technique, (3) improvement in 
the method of handling HIRS and MSU limb correction, and (4) improvement in 
the method of accounting for MSU surface emissivity. These modifications as 
well as the implementation of the simultaneous solution physical retrieval 
technique have measurably improved the quality of TOVS retrievals over the 
complex terrain of the Alps. 
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I. 


INTRODUCTION 


THIb grant was for partial support from the National Aeronautics and 
Space Administration to improve satellite vertical sounding retrievals over 
complex terrain conditions aB they related to ALPEX. Those features which are 
unique to complex land conditions are (1) the variable elevation of the 
terrain, (2) the large differences between surface air and surface skin 
conditions, and (3) the variability of surface emissivity at microwave wave- 
lengths. In the following section, each of these three points will be men- 
tioned very briefly. The reader will be referred to specific appendices for 
more detailed explanations. 

II. DISCUSSION 
A. SOFTWARE DEVELOPMENT 

A major development undertaken through this contract was the acquisition 
and application of high resolution topography on the McIDAS system. The 
previous 100 km resolution terrain file was improved to 15 km, which is better 
than the resolution of the HIRS instrument and therefore adequate for meso- 
scale sounding applications. The high resolution topography is necessary for 
evaluating the surface contribution term of the radiative transfer equation 
which is required for the proper interpretation of the radiance measurements 
in terms of atmospheric profiles. It is also important in establishing a 
reference level for the pressure height profile. 

A new technique for finding the surface skin temperature in the presence 
of cloud and reflected sunlight was implemented in the ALPEX retrieval soft- 
ware. The method uses the three window channels for 2 FOV and involves the 
minimization of a function relating surface temperature, cloud temperature, 
cloud amount and reflectance (Smith, et al., 1983; Appendix A). 


The statistical limb correction used in operational TOVS processing is 
known to be deficient. Because a random sample of cloud conditions are 
included in the regression dependent data base, the statistical correction 
tends to overestimate the limb correction for cloudy fields of view and 
underestimate it for clear fields of view. As a consequence, the statistical 
limb correction produces noise in the adjacent cloudy and clear fields of 
view. The ALPEX software was modified to account for the angular dependence 
of the radiances measured in the IR and MSU channels in a physical (rather 
than statistical) manner (Appendix A and C, Smith, et al., 1983 and LeMarshall 
and Schreiner, 1985). 

An algorithm for calculating the surface emissivity for the microwave 
channels has also been developed for the ALPEX algorithm. A detailed descrip- 
tion of the algorithm and evaluation is given in Appendix A. 

In all of the ALPEX processing it was intended that we maintain a 
standard benchmark version of the retrieval system against which to measure 
improvements. The benchmark version was essentially an operational iterative 
retrieval algorithm which does not Include the refinements carried out under 
this grant. Our four principal ALPEX orbits were processed with the benchmark 
version and reprocessed with several of the improvements. A systematic 
evaluation of each of the improvements was not performed, although comparisons 
of the benchmark version versus the "new and improved" iterative algorithm are 
described in Appendix A. 

In addition to the "new and improved" Iterative algorithm, the ALPEX 
grant was used in developing the "matrix inversion" retrieval algorithm 
(Smith, et al., 1984; Appendix D) . Briefly, this algorithm provides a simul- 
taneous solution for temperature and water vapor profiles and surface skin 
temperature from radiance observations. This approach has been shown to be 
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more Independent of the first guess than the iterative algorithm and also 
to yield greatly Improved moisture profiles. 

A feasibility study comparing single field of view (FOV) TOVS retrievals 
versus the box of 3x3 fields of view most TOVS software uses. The advantage 
of one FOV soundings would be in describing the atmosphere on a mesoscale 
level in the mountainopjs terrain. It waB found that the single FOV was most 
advantageous at the borders of mountains In two coses. The single FOV 
retrievals showed close agreement to 3x3 retrievals in areas of homogenous 
terrain. Results of this experiment are detailed in Appendix B (Loechner, 
1983). 

As part of the verification effort, software was developed to process 
ECMWF data (level II and III) on the McIDAS. These data can also be used to 
provide a first guess temperature profile though to date all retrievals have 
been made using a statistical first guess based on the NESDIS operational 
retrieval coefficients. 

It should be noted that the new ALPEX related algorithm development has 
been largely accomplished on the new McIDAS IBM 4381. The performance of this 
contract has provided the major impetus for transferring the TOVS man-inter- 
active software to the new machine and thereby is providing continuity in the 
SSEC satellite retrieval research effort. 

B . DATA 

Level 1-B data containing the raw radiances for all HIRS and MSU channels 
for NOAA-6 and 7 from the period 1 March - 30 April 1982, which corresponds to 
the first special observing period, were collected and preprocessed to produce 
earth located calibrated radiances. In addition, the ECMWF level Ill-a and 
level Il-a data for specific cases were received and are available on magnetic 
tape. The periods covered are: Level Ila; 1-10 March 1982 and 21-30 April 
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1982, Level Ilia; 1-9 March 1982 and 21-26 April 1982. In this case level Ila 
refers to a global observation data set for the area 30N-60N and 30W-37E. 

Level Ilia designates the operational global analysis for 0000, 0600, 1200 and 
1800 GMT on 1.875° latitude/longitude grid for the standard levels from 
1000 mb to 10 mb. 

In support of the research carried out under this contract, the University 
of Wisconsin sent out the following data sets to the various international 
ALPEX working groups. A total of eight NOAA-7 overpasses or four pairs were 
used. Two pairs, each covering the continent of Europe at the nominal time of 
4 March 1982 1300 GMT and 5 March 1982 0300 GMT; and two pair covering the 
Tasman Sea from the east coast of Australia to the International Date Line for 
the times 0300 GMT and 1500 GMT on 28 October 1982. In addition, software and 
coefficient files to read and process the ingested NESDIS TOVS l-B data sets 
were included. 

Meteosat data covering the period 4-5 March 1982 were requested of the 
European Space Agency and received. From these data cloud drift and water 
vapor winds were derived and are available on magnetic tape. 


III. CONCLUSIONS 

The purpose of this two year grant was to develop retrieval algorithms to 
deal with the unique problems of satellite retrievals over complex terrain. 

To a large extent this has been accomplished. In addition, this grant pro- 
vided the impetus for transferring software from the McIDAS/Harris system to 
the McIDAS/lBM-4341 and, most recently the IBM 4381 network. Also, Inter- 
national TOVS Processing Package (ITPP) versions I and II have been made 
available to users around the world. Finally, much of the ’•esearch work 
performed under this grant at the University of Wisconsin was in support of the 
First and Second International TOVS Conferences held in Igls, Austria. 
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Dota sets using both algorithms (1TPP versions I and II) have been gener- 
ated and statistics calculated comparing the retrieval results to ECMWF analy- 
ses. The results will be available in the Techulcal Proceedings o£ the Second 
International TOVS Conference. Cloud drift winds as well as water vapor winds 
for Meteosat data over the Alpine region nn March A and 5, 1982 have been 
generated and are available to Interested users. 


5 



REFERENCES 


LeMarshall, John F. and Anthony J. Schreiner, 1985: "Limb Effects in Satellite 

Temperature Sounding". Submitted to JCAM, 

Loechner, Franh 1983: "Temperature Retrievals over Mountains using one FOV". 

The Technical Proceedi ngs of t hi? First International TOVS Study 
Conference . Ed. Menzel, CIMSS. 132-141. 

Smith, W. L., H. M. Woolf, C. M. Hayden, A. J. Schreiner, J. F. LeMarshall, 
1983: "The Physical Retrieval TOVS Export Package". The Technical 

Proceedings of the First International TOVS Study Conference . Ed. Menzel, 
CIMSS. 227-278. 

Smith, W. L . , H. M. Woolf, and A. J. Schreiner, 1984: "Simultaneous Retrieval 

of Surface and Atmospheric Parameters - A Physical and Analytically Direct 
Approach". Presented at the Workshop on Remote Sensing Retrieval Methods. 
Williamsburg, VA, 30 Oct. - 2 Nov. 1984. 


4/AJS/04 


6 



Appendix Ai Physical Retrieval TOVS Export Package 
W. L. Smith, et al. 

excerpted from 

The Technical Proceedings of the First International 
TOVS Study Conference (editor P. Menzel) 



I 

\ 


THE PHYSICAL RETRIEVAL TOVS EXPORT PACKACE 

W.L. Smir.h*, H-'.'i. Woolf*, C.M. Hayden*, A.J. Schreiner^, J. F. LeMarahall^ 
Cooperative Institute for Meteorological Satellite Studies (CIMSS) 
University of Wisconsin 
Madison* Wisconsin 

1 NOAA/NESDIS, Development Laboratory 

2 Space Science and Engineering Center 

3 Visiting Scientist* Bureau of Meteorology, Melbourne, Australia 
1 • Introduction 

The satellite sounding radiance data received at the University of 
Wisconsin is converted into meteorological parameters using a physical 
solution of the radiative transfer equation, rather than empirical 
regression, in order to account for the influences of surface variables 
(i.e., terrain elevation, emissivity, and temperature) in the profile 
solutions and t.o permit the infrared data to be utilized for overcast as 
well as clear or partly cloudy sky conditions. These attributes of the 
physical solution are particularly important for the CIMSS applications of 
the data to real-time mesoscale weather analysis and forecasting for the 
North American region. 

The initial "TOVS Export Package" developed by the CIMSS, and made 
available to numerous direct readout data users, has been modified to 
Include the physical retrieval algorithm. As in the previous version, the 
new package produces temperature and moisture profiles from a 3 x 3 array 
of HIRS radiances and spatially Interpolated MSU observations. Thus, the 
resultant linear resolution and spacing is about 75 km. In addition, the 
new Physical Retrieval TOVS Export Package provides a measure of total 
ozone concentration based upon the stratospheric carbon dioxide and ozone 
channel (emission to space) radiance observations. Advantages of the phys- 
ical retrieval algorithm over the prior statistical regression method are: 

a. surface parameters such as skin temperature, emissivity, and elevation 
are treated as variables and therefore accounted for explicitly in the 
profile solutions; 

b. the reflected solar contribution to the shortwave window channels is 
determined and accounted for in the solution for surface temperature; 

c. the effect of clouds on the infrared channel observations is 
determined, thereby permitting the use of an optimum number of 
Infrared observations under all cloud conditions; 

d. the radiance data is not "limb-corrected", thereby eliminating a 
potential source of error in the data processing (the preprocessor 
(PROGRAM TOVPRE) now permits the user to specify whether or not limb 
corrections are to be made); 
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e. the solutions are in no way related to radiosonde data ao that their 
accuracy is ouch more Independent of geographical location and 
observation time than statistical regression solutions (with the 
exception of the regression first guess option). 

In this paper, prepared for the participants of the First International 
TOVS Study Conference, the physical algorithm is presented, its perfor»ps.ni:* 
on the Alpex and Tasman Sea NOAA-7 orbits chosen for the conference is 
described, and the documentation needed to implement the new software is 
given. 

2. Physical Retrieval System 

This document assumes the reader is already familiar with the observational 
characteristics of the TOVS; if not, the reader is referred to the articles 
provided in the bibliography. As a quick review, Table I shows the 
characteristics of the TOVS (H1RS and HSU) channels. Figure 1 shows a 
schematic of the TOVS Physical Retrieval Processing System. The ingest, 
calibration, and earth location are the some as given in the original 
version. In the physical retrieval system there is no need to angle or 
"limb-correct" the measurements since the weighting (i.e., atmospheric 
transmittance) functions can be computed for the angle of observation as 
described in section 2 i, 

a. First guess profiles 

j/nysical retrieval algorithm solves for a departure of the true 
atmospheric temperature, moisture, and ozone profiles from a "guess'* 
condition as a function of the departure of the observed radiances from 
those calculated for the "guess" profiles. The guess options included are: 
a latitude and time interpolation of the ten "Standard Atmospheres" (one 
for each of five latitude zones in winter and summer) , by regression 
estimates from the non-cloud-contaminated radiances (microwave and 
stratospheric infrared); estimates produced by a dynamical forecast model. 
The first two options are available in the physical retrieval export 
package; Bince accessing forecast-model output is highly machine- and 
installation-specific, each user must provide his own interface. 

b. Selection of "clearest" radiances 

The profiles are produced from 3x3 arrays of 26 kilometer resolution HIRS 
f ields-of-view (the courser resolution microwave data having been 
interpolated to the HIRS locations). The "clearest" observations are 
selected on the basis of the 11 pm window channel brightness temperature 
being within 3°K of the maximum observed amongst the 9 measurements. The 
average of all the "clearest" observations is generated for the sounding 
retrieval process. For the microwave channels and the stratospheric 
sounding infrared channels, the average is taken over all nine fields-of- 
view without regard to the 11 pm window brightness temperature. A check is 
made to determine whether the "clearest" observations are too clouu 
contaminated to afford an accurate profile retrieval. Using the 
11 pm brightness temperature and the "guess" temperature profile, if the 
former is colder than the guess 500 mb temperature no retrieval 1 b 
attempted. This check also safeguards against the use of precipitation 
attenuated microwave observations in the profile solutions. 
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c. Surface data 


The accuracy of the retrieved profiles can be affected significantly 
through the use of surface data. The surface parameters required for 
radiative transfer calculations are surface pressure, surface skin 
temperature, and surface emlssivlty. The "shelter" temperature and 
humidity are desired as constraints on the retrieved temperature and 
humidity profiles. 

(1) pressure 

Surface pressure Is computed hydrostatically using an estimated 1000 mb 
height, an estimated 1000 mb temperature and assumed lapse rate of S^C/km, 
and an estimate of the surface elevation for the radiance observation 
location. That is 
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Where g is gravity, R is the gas constant for air, and y the assumed lapse 
rate. The 1000 mb height and temperature are estimated from either an 
analysis of station observations reduced to che 1000 mb level, a forecast 
(aB in the case of first-guess generation, the user must provide the 
interface to such data), or climatology. The height of the observed 
surface is obtained from 10 nautical miles horizontal resolution topograph- 
ical data (this high-resolution topography is, at present, available for 
the Northern Hemisphere only: elsewhere, the previously supplied 60 

nautical mile resolution information must be used). The elevation for each 
HIRS field-of-viev is obtained through a f ield-of-view response weighting, 
and subsequently the elevation for each f ield-of-view is spatially averaged 
in the same manner sg the radiances themselves, as described in (2b). 

(ii) aii temperature and humidity 

The surface air temperature is estimated from an analysis of station 
observations reduced to 1000 mb by using the assumed lapse rates for the 
reduction to 1000 mb and the actual height of the observed topography. 
Surface humidity is estimated from an analysis of surface dewpoint 
depression. If station data is not used, then the surface air temperature 
and dewpoint are left as unknowns in the profile retrieval process. 

(ill) 6kin temperature 

The surface skin temperature is estimated from both longwave (11 pm) and 
shortwave (3.7 pm and 4.0 pm) window channels of the HIRS. During the day, 
the shortwave windows must be corrected for reflected sunlight 
contributions prior to the surface temperature retrieval. Assuming that 
the surface reflectivity is constant within the 3.7 - 4.0 pm region it can 
be shown (Smith, 1980) that 


R(3. 7 pm) - B(3.7po,T B ) 
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provided that the atmospheric transmittance is the same for both spectral 
channels. In (2), R is the observed radiance for H1RS window channel, B 
the Planck radiance, and T * s the effective radiating temperature of the 
sun, assumed to be 5800°K. u The only unknown in (2) is the effective 
radiating temperature of the earth-atmosphere for the 3.7 - 4.0 pm region, 
T_. Rearranging (2) gives 

IS 

f (T ) - KB(4 . 0pm,T )-B(3. 7pm,T ) - KR(4.0pm)~R(3. 7pm) (3) 

D D D 

Thus, f(T R ) is calculated once and for all, and tabulated to permit rapid 
determination of T„ from actual observations of R(4.0 pm) and R(3.7 pm) 
observations (Channels 18 and 19 of Table 1). 

The effective surface temperature for the 4.0 pm and 11 pm regions is then 
calculated from the observed (solar corrected) brightness temperatures 
using the relation 
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where T fi refers to the^brightness temperature computed for a "guess" 
surfacetemperature, Tg , and atmospheric temperature and moisture profile 
condition, and t(p ) is the transmittance of the atmosphere above the 
surface pressure level, p . Equation (4) is derived from a form of the 
radiative transfer equation in which the Planck radiance dependence upon 
temperature is linearized using a first order Taylor expansion. Since (4) 
is not exact, the solution is Iterated several times until | T n+ - T | 

< 0.01°K. The initial value for T is obtained from the "split window® 
approximation, T ' * « Tg (11 pm) + S 0 . 5 [ T_ (11 pm) - T fi (8 pm)). As shown in 
Figure 1 , the temperature and moisture profiles usedior the final 
determination of T are provided by the microwave brightness temperatures 
and H1RS water vapor channel radiances as described later. The longwave 
(11 pm) and shortwave (4 pm) surface temperature values are checked for 
consistency and used to diagnose the presence of cloud contamination. 
Provided the results are reasonable for cloudfree sky conditions, they are 
used separately for the calculation of longwave (15 pm) and shortwave 
(4.3 pm) radiances in the infrared temperature profile retrieval process. 
The differences between the two surface temperatures are attributed to 
surface emissivity differences for the two spectral regions. 
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(lv) microwave emissivity 


Unlike the emissivity of the earth' a surface at infrared wavelengths, the 
surface emisslvity at microwave frequencies can depart drastically from 
unity, depending upon surface moisture conditions. Thus, the surface 
emissivity must be determined and used explicitly in the calculation of 
outgoing radiance. The microwave "window" channel at 50 GHz (Table 1) was 
included in the TOVS specifically for determining the effects of surface 
emisslvity and liquid water on the microwave sounding channel radiance . 

It can be shown that the surface emissivity, c^, is given by the expression 
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where t(p) is the transmittance of the atmosphere above pressure level p. 
For the microwave window region (50 GHz), the atmospheric (integral) 
quantities are small compared to the observed brightness temperature and 
surface term; thus the result is weakly dependent upon the temperature 
profile condition assumed. If a poor first guess temperature profile 
condition is used (e.g., climatology) the entire surface emissivity and 
microwave temperature profile retrieval process i6 repeated. 


d. Microwave temperature profile retrieval 

The temperature profile is calculated from the microwave brightness 
temperatures observed in MSU channels 2,3, and A, using the iterative 
solution (Smith, 1970) 

I 1W - T b" (V 1 )IW< V p) 

t”*‘<P> - l"(p) + — J 

l W(v p) 

i«2 x 


where Tg” is the brightness temperature calculated from the temperature 

profile T n (p), and W(v 1 ,p)“dx(v 1 ,p)/dp, In the calculation of T n (p), the 

surface emisslvity calculated from the 50 GHz brightness temperature 
measurement, the surface temperature given by the H1RS 8 vm and 11 ;«n 
brightness temperatures using the "split window" approximation, and the 
surface pressure obtained from the high resolution topographical data are 
used as boundary conditions. The solution for T(p) is iterated until the 
condition 

l T B n<V i > " T B n ” l<V i ) I < °* 05 ° K 
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1b satisfied for all channels. If surface air temperature observations are 
being used, the profile below the 700 mb pressure level is adjusted in a 
linear fashion to blend into the observed air temperature at the surface. 


e. Initial water vapor profile estimate 


The initial profile of water vapor is derived using the mean relative 
humidity solution given by Smith and Zhou (1980). Assuming the relative 
humidity is constant within each layer sampled by the H1RS water vapor 
channels (10, 11, and 12), then 
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where p(v.) is the pressure level of maximum water vapor sensitivity for 
the watervapor channel v., . Rh° is a guess layer relative humidity value (a 
value of 0.5 is used initially), R°(v.) is the radiance calculated for the 
guess humidity condition, and W B (v^) is the water vapor weighting function 
given by 


W°(v 1 ,p) * g 


U°(p) 8t(v 1 , P ) 3B(v^,T) 

q°(p) &p 3p 


where U°(p) is the integrated precipltable water above the pressure p, and 
q e (p) is the mixing ratio as implied by the guess relative humidity 
condition, and the atmospheric temperature profile (which, at this point, 
has been derived from the MSU brightness temperature observations). For 
the HIRS water vapor channels 10, 11, and 12, the pressures of maximum 
sensitivity are assumed to be 800, 600, and 300 mb, respectively. Given 
these three relative humidity values from (7) a profile of relative 
humidity is constructed by linear interpolation with respect to pressure. 

If surface data Is available, it is used to Interpolate between the surface 
and the 800 mb level, otherwise the relative humidity is assumed to be 
constant below the 800 mb level. The relative humidity above 300 mb is 
also assumed to be constant. The solution given by (7) 1 b repeated after 
the infrared brightness temperatures are used to derive an improved 
estimate of the temperature profile, using the relative humidity values 
obtained initially. Although the layer mean relative humidity 
approximation is not an exact simultaneous solution for the relative 
humidity profile from the three water vapor channel radiances, it is a 
stable approximation in that the influence of cloud or measurement errors 
is minimal. 

It should be noted that, even if the atmosphere is cloudy, a relative 
humidity profile is derived. In this case, the solution is produced by 
using the microwave temperature profile and the 11 pm brightness tempera- 
ture to estimate the cloud pressure and temperature. These cloud 
parameters are then used in the calculation of R° (v ). The relative 
humidity results below the cloud level are a consequence of linear 
interpolation between the relative humidity sensed above the cloud level 
and the observed surface value. 
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f. Infrared temperature profile retrieval 

After deriving the water vapor profile, the surface skin temperatures for 
the short and long wavelength temperature sounding channels are obtained by 
iterating (4). These together with the microwave specified temperature 
profile are used to calculate the "expected" values of the brightness 
temperatures observed in the 4.3 pm and 15 pm temperature sounding channels 
(Table 1). If the visible channel reflectance is less than 25 % and the 
longwave surface temperature is greater than a threshold of flv? degrees 
below the observed or microwave retrieved surface air temperature, then all 
the infrared channels are presumed to be free of cloud contamination, and 
are used in the temperature profile retrieval process. If the "clear test" 
fails, then an attempt 1 b made to determine which of the infrared channels 
are free of cloud contamination. By assuming that cloud6 always reduce the 
brightness temperatures by a magnitude dependent upon the transmissivity of 
the spectral channel, all channels of equal or greater opacity than a given 
channil are presumed to be clear if the deviation between the observed and 
calculated brightness temperature for that channel 16 greater than or equal 
to zero. If HIRS channel 7 satisfies this criterion, all shortwave and 
longwave channels are presumed to be unaffected by the cloud and used in 
the retrieval process. If, however, this criterion 1 b not satisfied by 
channel 7, then only the longwave channels satisfying the criterion are 
used; the shortwave channels are not used because of the difficulty in 
accurately correcting these measurements for cloud-reflected sunlight 
contributions. 

Having determined the infrared channels unaffected by cloud contamination, 
the temperature profile is given by the iterative solution 

l [T B ( Vl ) - T B n (v i )]W(v i ,p) 

T n+1 (p) - T n (p) + ^ (8) 

“ WCv^P) 

i“l 


where m is the number of cloud free infrared channels and the weighting 
function is given by the relative contribution function 


W(v i p) - 
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If^any longwave channel, v^, is presumed to be cloud contaminated, then 
Tg (v^) is presumed to be equal to Tg(v.) in equation (8), in order to 
blend 1 the solution into the microwave retrieval below the cloud level. The 
iteration is terminated as soon as the criterion 


t b D< v 1 ) - T B n ~ 1<v i ) < °* 05 ° K 


233 



Is satisfied for all channels. The temperature profile retrieved from the 
infrared brightness temperatures, using the microwave result as the Initial 
guesB, Is the "final" estimate. 

g. Final water vapor profile estimate 

The final water vapor profile ia calculated using (7) with the infrared 
temperature profile and the final surface temperature estimates being used 
in the calculation of R°(v ) and W°(v,p). The Rh® values ore taken as the 
estimates initially derives on the basis of the microwave temperature 
profile. 

h. Ozone retrieval 

The total ozone concentration is retrieved from the ozone channel 
brightness temperature using the method of Ma, Smith, and Woolf (1984). In 
a manner analogous to the water vapor solution, the ozone concentration is 
given by the inverse solution of the radiative transfer equation. 
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°W°(v 1 ,p)dp 


where is channel 9 (Table I) and 

0°(p) 8t°(v p> 3B(v T) 

W (v ,p) - g — * - 

m Ip) dp dp 

0(p) is the ozone concentration profile, and m(p) is the ozone mixing ratio 
profile. The initial ozo!.«e mixing ratio and concentration profile is 
obtained using regression relations developed from a large sample of 
ozonesondes and radiances for the stratospheric H1R5 channels (1-4) 
synthesized by radiative transfer calculations, from the ozonesonde 
temperature profile data. Thus the temperature profiling channels are used 
to specify the gross vertical structure of the ozone profile (via 
temperature-ozone correlation), and the ozone channel radiance is used to 
estimate the absolute magnitude of the total concentration. If a cloudy 
condition prevails, the surface pressure used in (9) is the cloud pressure 
determined from the 11 pm window brightness temperature and the final 
temperature profile. 

1. Calculation of transmlttances 

An essential ingredient of any physical retrieval scheme i3 the 
availability of accurate and computationally efficient ("fast") procedures 
for calculating profiles of atmospheric transmittance in the various 
spectral intervals represented by the infrared (HIRS) and microwave (MSU) 
sounding channels. These procedures must account for variations in 
spectral response, atmospheric composition and state (temperature and 
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moisture dependence)* and departure from vertical viewing (slant path or 
zenith-angle dependence). 

The requirement for computational speed precludes the employment* In 
routine processing, of so-called "exact” (in the spectroscopic sense), or 
line-by-line, methods. Therefore, models have been developed (McMillin and 
Fleming, 1976; Fleming and McMillin, 1977; Welnreb and Neuendorffer , 1973; 
W'ilnreb et al. , 1981) to satisfy the dual requirements of speed and 
accuracy — the latter being defined in terms of degree of agreement with 
exact calculations. The tirst three references cited above describe the 
algorithms developed for the infrared; the first relates these specifically 
to HIRS, and also discusses the method U6ed with MSU. 

When the design decision was made to apply the physical retrieval algorithm 
to non-limb-corrected TOVS data, it was determined that the existing MSU 
transmittance model would not suffice. For an arbitrary scan angle, it 
entails calculation of transmittance profiles for a number of angles, 
followed by convolution with the antenna gain pattern — a computationally 
costly process. Therefore, a variation of the model was developed, in 
which a separate set of regression coefficients, relating absorption to 
temperature and moisture, is generated for each of six equally spaced scan 
angles. The influence of the antenna pattern is incorporated into the 
absorption profiles used for the regression. Coefficients for arbitrary 
scan angles are obtained dynamically by linear interpolation. The model 
has been verified by comparing antenna temperatures calculated from "exact" 
and regression-estimated transmittance profiles; agreement is very good. 

Transmlttances for the ozone retrieval are computed using the method of 
Gruenzel (1978), with contributions from water vapor lines and continuum 
obtained by the procedures described in Weinreb, et al. (1981). 

The transmlttances are adjusted in order to minimize differences between 
observed and calculated radiances. The adjustments are made assuming that 

lnT^p) ■» y(v i )lnt(v i ,p) 

where yCv.^) is the adjustment factor and t(v.,p) is the unadjusted 
transmittance. The adjustment factors have Been determined on the basis of 
a very large sample (several thousand) of "coincident" radiosonde and 
radiance observations. The y*s are calculated so as to minimize the RMS 
difference between observed and calculated radiances for all atmospheric 
conditions ar.d latitude zones. The values assumed for all NOAA satellites 
in the TOVS Export Package are given in Table 2. 

3, Results From NOAA-7 ALPEX Orbits 


The Physical Retrieval Export Package was applied to the NOAA-7 orbits 
chosen for the TOVS Study Conference. The program options used were 
"Regression Guess" and "k.'th Surface Data". The regression coefficients 
used were the NESD1S operational for the overcast cloud (Path 3) condition. 
The surface data was obtained from station reports for the ALPEX orbits and 
from Australian Bureau of Meteorology analyses for the Southern Hemisphere 
orbits. The effect of the different program options upon the results are 
discussed in section 5 of this report. 
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The actual retrieval values and radiosonde data for the Alpex orbits are 
given in an Appendix. Also, the retrieval values and contour analyses for 
the Southern Hemisphere orbits chosen for the Study Conference tire given in 
an Appendix, The geopotential thickness values are for the layers 
specified in the Study Conference invitation and ore given in decameters. 

In this Bectlon we present comparisons of the various TOVS geopotential 
thickness analyses with those of radiosonde data produced by the European 
Center for Medium range Weather Forecasting (ECMWF) for the ALPEX orbits. 
The analysis of the TOVS retrievals is produced by Barnes interpolation 
(Barnes. 1973) at a spatial resolution of 75 Km. The ECMWF analyses are 
based on optimal interpolation (Bengtsson, 19P1) at a spatial resolution of 
200 Km. The different resolutions are consistent with the different 
horizontal spacings between TOVS and Radiosonde sounding data (see the 
Appendix) . 

Figures 2-6 show the analysis comparisons for gcopotential thickness. As 
may be seen from these figures, the patterns of the two different analyses 
produced from TOVS (SAT) and from Radiosonde (ECMWF) data are nearly 
Identical, however, there are significant differences in the horizontal 
scale and gradients resolved by the two data/analysis systems, In 
particular, smaller horizontal scales seem to be resolved in the TOVS 
analyses. For example, the shortwave ridge in the 1000-700 mb thickness 
along the east-coast of Spain in the 03Z TOVS analysis (Figure 2c). Also, 
the amplitude of the trough and ridge axes over the Oceanic areas are 
definitely greater in the TOVS analyses; for example, the ridge axis over 
the ocean between Spain and England in the 1000-500 mb thickness at 03Z on 
5 March (Figure 3c). On the other hand, at upper tropospheric levels there 
is a definite weakening of the horizontal gradient of TOVS derived 
thicknesses relative to the radiosonde data. For example, the horizontal 
gradient of the 300-100 mb thickness across Italy for 00Z on 5 March is 
about 180 meters for the ECMWF analysis, compared to about 90 meters for 
the SAT analysis. The magnitude of the high over the Mediterranean also 
differs by about A0 meters, the SAT analysis being weaker than the ECMWF 
analysis. The loss of horizontal resolution of the TOVS soundings relative 
to radiosondes with altitude is due to the loss of vertical resolution of 
satellite soundings with Increasing altitude. The relatively low vertical 
resolution of the satellite soundings in the tropopause region leads to 
large bias and relative differences with radiosonde observations (section 
A). 

Figure 7 shows the twelve hour change of precipitable water above 1000, 

700, and 500 mb levels as retrieved from the TOVS data. One can see 
good continuity with time and with the thickness patterns displayed in 
Figures 2-6. 

Finally, Figure S shows the total ozone retrieved from the N0AA-7 TOVS 
observations for A and 5 March, 1982. Unfortunately, the ground truth data 
was not available to verify the absolute accuracy of these retrievals. Of 
particular interest is the very rapid movement of the high ozone center 
over the English channel at 13Z on A March to the Mediterranean at 03Z on 5 
March. This movement, and in fact the entire ozone pattern, is highly 
correlated with the patterns and time displacements seen in the 300-100 mb 
thickness shown in Figure 6. This correlation is presumably due to the 
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relation between ozone concentration and geopotential thickness In this 
layer to tropopause height and temperature. High stratospheric ozone 
concentrations should be associated with low tropopause height and low 
tropopauae height are associated with warm tropopause temperatures or 
geopotential thicknesses. Thus, there is a strong positive correlation 
between the geopotential thickness of the tropopause layer (l.e., 300-100 
mb) and the stratospheric ozone concentration. This correlation 16 clearly 
seen in the TOVS retrievals (compare Figures 8 and 6). 

A. Error Statistics 

Figure 9 presents bias and standard deviations between the analyses of TOVS 
retrieved temperatures and radiosonde temperature observations (ECMWF) as 
well as the bias and standard deviations of each of these analyses with 
respect to actual radiosonde observations. Because of occasional large 
digital transmission errors in the radiosonde data, observations which 
exceeded 2® (Figure 9c) or 5° (Figure 9d) deviation from the ECMWF analysis 
were discarded from the error statistics. No filtering of TOVS data was 
performed. The reader is reminded that the TOVS retrievals were produced 
using a stratospheric H1RS and MSU channel regression "first guess" and 
incorporated surface data. The TOVS and ECMWF analyses were performed with 
grid spaclngs of 75 km and 200 km, respectively. 

As may be seen by comparing Figure 9(b) with 9(c) or 9(d), the bias and 
standard deviations between radiosonde observations and the TOVS and ECMWF 
analyses are comparable, except in the tropopause regions. This result is 
somewhat surprising since the TOVS analyses are completely Independent of 
the radiosonde data whereas the ECMWF analyses are completely dominated by 
radiosonde observations over the European region studied. One explanation 
for the comparability of the magnitudes of the standard errors is (a) the 
lack of small horizontal scale detail in the ECMWF analyses (i.e., 
horizontal smoothing) and (b) the inherent vertical smoothing of the TOVS 
profiles relative to radlsonde observations. The much larger bias error in 
the TOVS retrievals in the tropopause region merely demonstrates that 
vertical smoothing is more detrimental than horizontal smoothing in this 
region of sharp vertical temperature variation. The combination of the 
error statistics Bhovn in figure 9 indicate that the TOVS level temperature 
retrieval errors are generally less than 2®C. Errors in thickness 
temperatures should be less than the error of level temperature values 
because of the vertical smoothing inherent in the TOVS retrievals. 

5. Impact of Profile Guess and Surface Data 

The influence of profile first guess and surface data options of the 
retrieval program were studied using the western orbit of A March, 1982. 
Figure 10 shows the bias and standard deviations of the TOVS level 
temperature values with respect to the ECMWF analyses for four different 
program options; (a) regression guess with surface data, (b) climate guess 
with surface data, (c) regression guess without surface data, and (d) 
climate guess without surface data. As may be seen, the main influence of 
the various options is upon the bias error distribution of the retrievals. 
The major influence of surface data is below the 700 mb level where 
standard deviation and bias error Improvement of more than 1°C can be seen. 
The influence of the profile first guess (regression vs. climate) seems to 
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greatly Impact the bias error for the case where surface data Is utilized. 
This may be due to the structure differences created by blending the "first 
guess" profile into the surface observations. Surprising is the fact that 
without surface data (Figures 10c and lOd), little difference is produced 
by the profile first guess. Comparisons of horizontal analyses of the 
level temperatures for the various first guess options (not shown) revealed 
little difference in the patterns but significant differences in the 
absolute temperature values as suggested in Figure 10. The authors believe 
that the most reliable retrievals are produced for the regression gues6 
with surface data option. Certainly, surface observations are extremely 
valuable when dealing with complex terrain conditions. 

6. Summary and Future Improvements 

The physical retrieval "Export Package" ha6 been described. Analyses of 
its performance on the NOAA-7 orbits chosen for the TOVS study conference 
Indicate satisfactory performance. A major Improvement yet to be 
implemented is the use of full resolution AVHRR data to better define 
surface skin temperature and to account for the influence of clouds upon 
the TOVS sounding channel radiances. An effort is now underway at the 
C1MSS to implement the AVHRR option to the TOVS sounding processing system. 
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TABLE 1. Characteristics of TOV Sounding Channels 
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Table 2. Transmittance - adjustment exponents ("gammas") 
for use in SUBROUTINE TAUGAM 


H1RS MSU 


Channel 

Gamma 

Channel 

Gamrw 

1 

1.069 

1 

1.000 

2 

0.942 

2 

1.019 

3 

1.108 

3 

0.999 

4 

1.031 

4 

0.899 

5 

0.972 



6 

0.965 



7 

0.909 



8 

1.000 



9 

1.000 



10 

1.120 



11 

1.032 



12 

1.035 



13 

1.115 

• 


14 

1.028 



15 

1.057 



16 

1.002 



17 

1.480 



18 

1.000 



19 

1.000 
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Figure 1. Schematic of the TOVS Physical Retrieval Processing System. 
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Figure 6. Analyses of TOVS (SAT) and Radiosonde (ECMWF) 300—100 mb geopotential thickness (decameters) 




(a) WI. WHrtd «so*c tooonsicmioi h nw UK |» Mi 


•HI 



me. MAIM H80VC SaOHfli:r*lOO) A NRR JBBJ 13J MI (f) me. MMCR AB0IC 9«£tBtCflil00> 9 HAR 1887 037 MI 


Figure 7. Analyses of TOVS (SAT) Total Precipitable Water for the Atmospheric Column 
above (a) and (b) 1000 mb, (c) and (d) 700 mb, and (e) and (f) 500 mb. 
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Figure 10, Bias and Standard Deviations between ECI'tWF analyses of radio- 
sonde data and TOVS level temperature values produced for 
four different program options, (a) regression guess with 
surface data, (b) climate guess with surface data, (c) regression 
guess without surface data, (d) climate guess without surface 










Appendix to Smith et 

The following pages contain plots of TOVS and radiosonde geopotential 
thickness ( decameters ) and precipitable water ( cm x .100 ) for the 
TOVS Study Conference N0AA-7 orbits. For the southern hemisphere, TOVS 
data and analyses are also given. Documentation for the TOVS Export Package 
software and data are included at the end of this Appendix. 
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DATA FROM THE TIROS-N/NOAA SERIES OF POLAR-ORBITING 
THE SOFTWARE AND SUPPORTING DATA CONTAINED IN SUBSE 
TAPE. THE ALGORITHMS AND DATA-PROCESS ING TECHNIQUE 
DEVELOPED FOR ’KC1DAS* (MAN-COMPUTER INTERACTIVE DA 
THE SPACE SCIENCE AND ENGINEERING CENTER OF THE UNI 
IN MADISON (UU/SSEC). INITIALLY BASED ON THE HARRI 
COMPUTER* MCIOAS HAS RECENTLY BEEN IMPLEMENTED ON T 
SOFTWARE HEREIN SUPPLIED HAS BEEN CONFIGURED FOR AN 
PROVIDE THE MOST PORTABLE SYSTEM FEASIBLE* GIVEN TH 
TO THOSE RESPONSIBLE FOR MAINTENANCE OF THE SYSTEM. 


PROCEDURES FOR ESTAB- 
CHIRS*MSU) SOUNDING 
SPACECRAFT* USING 
OUENT FILES ON THIS 
S WERE ORIGINALLY 
TA ACCESS SYSTEM) AT 
VERSITY OF WISCONSIN 
S ’SLASH-6* MINICOM- 
HE IBM 4341. THE 
IB M -OS SYSTEM TO 
E RESOURCES AVAILABLE 


ftid' 

QUESTIONS AND/OR COMMENT S CONCERNING THE SYSTEM AND ITS IMPLEMENTA- 
TION* ESPECIALLY IN RCGARO TO POSSIBLE ERRORS. SHOULD BE ADDRESSED TO 
MR. HAROLD M. WOOLF 

NOAA-NESDIS DEVELOPMENT LABORATORY 
UU/SSEC* ROOM 219 
1225 WEST DAYTON STREET 
MADISON. WISCONSIN 53706 

TELEPHONE ... 1608) 264-5325 
. <910) 266-2771 
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DNEXPORT 




NOTES ON FORTRAN SOURCE (FILES 2 AND 3> 

< • 1 1 f 

THE FIRST CARD OF EACH MAIN PROGRAM IN FILE 2 IS A FORTRAN COMMENT OF 
THE FORM 

C MAIN PROGRAM *XXXXXX* 

EACH MAIN PROGRAM* SUBPROGRAM « OR COLLECTION OF RELATED SUBPROGRAMS* 
CONTAINS A COMMENT CARD NEAR ITS BEGINNING OF THE FORM 
C -*» VERSION OF OD KHM YY (DAY MONTH YEAR) 

USERS FOR WHOM THIS IS NOT THE FIRST PACKAGE SUPPLIED SHOULD FIND 
THESE ’VERSION-DATES* A CONVENIENT MEANS OF DETERMINING WHICH ROUTINES 
HAVE EEEN CHANGED SINCE THEIR INITIAL IMPLEMENTATION. ALTHOUGH IN THE 
PAST UPDATES TO THE PACKAGE HAVE BEEN PROVIDED IN THE FORM OF SELECTED 
SOFTWARE AND COEFFICIENTS* THE MOST RECENT CHANGES ARE TOO EXTENSIVE* 

AND THE NUMBER CF USERS WORLD-WIDE TOO GREAT* TO PERMIT SUCH ’LIMITEO- 
EC1TION* UPDATING. HEREAFTER* UPDATES - TO SOFWARE* COEFFICIENTS* OR 
BOTH - WILL BE IN THE FORM OF A NEW TOTAL-PACKAGE TAPE. IT WILL BE THE 
RESPONSIBILITY OF THE INDIVIDUAL USER TO DETERMINE WHICH ’NEW ITEMS* ARE 
RELEVANT TO HIS INSTALLATION. 

NOTES ON SUPPORTING-DATA FILES 

i * A i * 

THE COEFFICIENTS IN FILES 7 THROUGH 12 ARE DETERMINED INITIALLY FOR 
EACH SATELLITE* AND CAN BE CONSIDERED FIXED FOR ITS OPERATIONAL LIFETIME 
UNLESS MAJOR CHANGES OCCUR IN THE PERFORMANCE OF ONE OR MORE SPECTRAL 
CHANNELS . * THIS IS ALSO TRUE FOR THE COEFFICIENTS IN FILES 15 AND 16. 
THE RETRIEVAL COEFFICIENTS IN FILES 14 AND 17 ARE UPDATED WEEKLY BY 
NE SO I S OPERATIONS AND CAPTURED FOR INCLUSION IN THIS PACKAGE. ORBITAL 
ELEMENTS SUPPLIED IN FILES 5 AND 6 ARE ALSO UPDATED WEEKLY. THESE ARE 
NEEDED FOR NAVIGATION (EARTH-LOCATION) OF REAL-TIME DATA OBTAINED VIA 
DIRECT READOUT OF THE SPACECRAFT (BEACON). THE REAL-DATA COEFFICIENTS 
AND ORBITAL ELEMENTS ARE PROVIDED TO ASSIST IN THE INITIAL IMPLEMENTA- 
TION? THE USER MUST MAKE HIS OWN ARRANGEMENTS FOR CONTINUED ACQUISITION 
OF SUCH INFORMATION FOR REAL-TIME APPLICATIONS. 

OCCASIONALLY* SOME OF THE M.SU CALIBRATION QUALITY-CONTROL LIMITS IN 
FILES 7 AND 8 MAY NEED TO BE ADJUSTED TO AVOID DISCARDING GOOD DATA 
AS A RESULT OF CHANGES IN GAIN OR OTHER INSTRUMENT CHARACTERISTICS. 
*■«»:* 

CONSTRUCTING THE TOVS 0 AT A-PROCESSING SYSTEM 

**«:« 

THE SYSTEM MAKES EXTENSIVE USE OF DISK FILES* BOTH SEQUENTIAL AND 
DIRECT-ACCESS, FOR EFFICIENCY IN I/O OPERATIONS. TO ASSIST IN ESTAB- 
LISHING THE PERMANENT DATA FILES (ORBITAL ELEMENTS, COEFFICIENTS* AND 
TOPOGRAPHY), PROGRAMS HAVE BEEN PROVIDED (IN FILE 2) TO READ THE DATA 
F R 0 v TAPE AND WRITE TO DISK. 


PROGRAM 

• t f 

PUTELE 


TAPE FILE ( S > 

4 « f I'i ♦ ■* * 1 * V 

5*6 


DISK LUN 

•• * * r * * f » 
11*12 


TYPE-* 

li 9 * * 

D*D 


♦ ■**•*# V 
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26 
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NOTTS : 

l UN r FORTRAN LOGICAL UNIT NUMBER 
* S = SEQUENTIAL* D = DIRECT-ACCESS (IBM BDAM ) 

V PROGRAM *RTVLCFTP» IS PROVIDED TO PERMIT EXTRACTION OF RETRIEVAL 
COEFFICIENTS FROM THE GENERAL COEFFICIENT-DATABASE TAPE THAT IS AVAIL- 
ABLE FPOM NESDIS OPERATIONS UPON SPECIAL REQUEST. 

a** PROGRAM • RTVLCXTP* IS ANALOGOUS TO *RTVLCFTP» , BUT PROVIDES THE 
RECORD STRUCTURE REQUIRED FOR INPUT TO PROGRAM *FXTIRO*. 

♦ PROGRAM * SUBTOP * IS PROVIDED TO CREATE A SUBSET OF THE GLOBAL HJGH- 
RESCLUTION TOPOGRAPHY FOR USE IN PROGRAM *FXTIRO*. IT IS ASSUMED 
THAT MOST USERS WOULD NOT BE ABLE TO MAINTAIN THE ENTIRE DATASET ON 
DISK. SINCE IT CONTAINS 16200 RECORDS. EACH 2 68 BYTES IN LENGTH. 

TWO PROGRAMS. *RAOBHlRS« AND *RAOBMSU* (AND ASSOCIATED SUBPROGRAMS) 
ARE SUPPLIED TO DEMONSTRATE THE PROCEDURES FOR CALCULATING, FROM RADIO- 
SONDE TEMPERATURE AND HUMIDITY PROFILES. HIRS RADIANCES AND BRIGHTNESS 
(EQUIVALENT-BLACKBCCY) TEMPERATURES. AND *SU ANTENNA OR BRIGHTNESS TEMP- 
ERATURES, RESPECTIVELY. THE REMAINDER OF THE SOFTWARE CONTAINED IN 
FILES 2. 3, AND A IS FOR PROCESSING HIRS AND MSU DATA TO OBTAIN PROFILES 
CF ATMOSPHERIC TEMPERATURE. HUMIDITY, GEOPOTENTIAL HEIGHT, AND GtOSTRC- 
PHIC WIND, AND FOR DISPLAYING AND MANIPULATING THOSE PROFILES IN VARIOUS 
WAYS. ESSENTIAL TASKS APE * INGEST* , *PRE PROCESS I NG * , AND •RETRIEVAL*. 
THE TERM *RTC* MEANS * R AD I AT I VE -TR ANSFER COEFFICIENTS*. AND *LCC* STANDS 
FOR *LIMB-C0RRECTI0N COEFFICIENTS*. 

I? * V \£ \ 


* * V T + V 9 

INGEST 


V * 3 T * 

function: produce CALIBRATED, EART 
METRIC MEASUREMENTS FROM TOVS *TIP* D 
THE ESSENTIAL DIFFERENCE IS IN THE TY 
TO HANDLE - EITHER ARCHIVAL ('LEVEL 1 

LEVEL 1-E: DATA THAT HAS BEEN THRO 
NESDIS OPERATIONS* TOVS GROUND SYSTEM 
PUTER TAPE TO USERS. UPON REQUEST, BY 
SUCH TAPE CONTAINS. IN ADD 
IG RADIOMETRIC MEASUREMENTS 
i PARAMETERS RECUIRED TO TR 
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OF NESDIS 
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REPRESENT 


CAL IBRATI 


H-LOCATED HIRS AND MSU RADIO- 
ATA. TWO VERSIONS ARE PROVIDED? 

FE OF INPUT OATA THEY ARE DESIGNED 
-B * > OR DIRECT-READOUT (REALTIME). 

UGH PRELIMINARY PROCESSING BY 
, AND IS PROVIDED ON STANDAP D COM- 
THE SATELLITE DATA SERVICES BRANC 
ITION TO VALUES (IN DIGITAL COUNTS 
, E A RTH-LCCATI ON INFORMATION AND 
ANSFOP.M THE RAW DATA VALUES INTO 
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197 RADIANCE OR BRIGHTNESS TEMPERATURE. THREE PROGRAMS ARE SUPPLIED FOR 
we THIS TYPE OF OPERATION: 
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PROGRAM 

INPUT CS) 

U1 

.ftt** 

.«*«*, vi t « 

Ip 2 

INVTAP 

1-B TAPE CIO) 
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DATA CARD 

164 

TOVT AP 

1-B TAPE(IO) 
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''•DATA CARD 


OUT PUT t S ) < ) = LUN 

PP INTOUT OF DATA TYPE* TIMES. AND 
LOCATIONS FOR EACH FILE 

SELECTED HIRSfll) AND HSU C 12) DATA ON 
DISK} PRINTOUT OF RELEVANT INFORMATION 


167 

Tp 8 


170 


TOVING 


HIRS DISKIll) CALIBRATED* LOCATED DATA ON DISM20); 
MSU Cl SK ( 12) PRINTOUT 

HIRS R TC ( 1 3 > 
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S ‘o 9 
190 


t-^DATA CARD 


notes: 

' NUMBER OF FILES TO READS PROGRAM TERMINATES ON DOUBLE EOF IF ENC CUNT 
ERED BEFORE COUNT IS SATISFIED. 

PHYSICAL FILE NUMBERS OF HIRS AND MSU DATA* PLUS BEGINNING AND END- 
ING TIMES. 

*. SATELLITE NUMBER <1 FOR NOAA-7* 2 FOR NOAA-P)S REQUIRED BECAUSE 
INFORMATION IN TAPE HEADER RECORDS MAY BE AMBIGUOUS. 

DIRECT-READOUT: DATA OETAINED ON-SITE BY DIRECT DOWNLINK FROM THE 
SPACECRAFT. THE SSEC SYSTEM USES DATA FROM. THE VHF C137MHZ) BEACON. 
SUCH DATA CAN ALSO BE ACQUIRED BY EXTRACTING *TIP» FROM AN HR»T DATA- 
STREAMS CREATION OF THE NECESSARY SOFTWARE IS THE RESPONSIBILITY OF THE 
USER. PROCESSING OF C I R ECT -R E ADOUT DATA IS MUCH MORE COMPLEX THAN THAT 
OF LEVEL 1-B » SINCE THC USER MUST DO EVERYTHING - DECOMMUTATION OF RAW 
DATA* NAVIGATION OR EARTH-LOCATION* AND IN-FLIGHT CALIBRATION. TWO PRO 
GRAMS ARE PROVIDED: 

PROGRAM INPUT(S) OUTPUTCS) 
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PREING 

DATA T A PE < 9 ) 
DATA CARD 

DISK FILE CIO) CONTAINING 
HIRS AND MSU DATA 
PRINTOUT 

DECOHMUTATED 

INGTOV 

DECOM DATA(IO) 

• * ORE-fLEMdl *12 ) 
HIRS RT C (13) 
ING-PARAMC15) 

CALIBRATED, 

PRINTOUT 

LOCATED DATA 

ON DISKI20) 


2 00 

2j]jl NOTES: 

21.1; 2 FILE TO PROCESS* IF MORE THAN ONE ON TAPE* FLAG FOR DETAILED (OIAG 


203 NOSTIC) MSU PRINT. 

2|M * UPDATED BY PROGRAM • PUTELE * WITH DATA OBTAINED FROM DIRECT READOUT 

ajljs . USER SERVICES • 

2 06 

2A7 THE DATA IN FILE 20 WILL HAVE THE SAME FORMAT. REGARDLESS OF THE 

2|je SOURCE AND TYPE OF INGEST, THIS FILE SERVES AS INPUT TO THE NEXT STEP. 
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PREPROCESSOR 

» • * f. • + ** 1 

FUNCTION: TRANSFORM CALIBRATED* EARTH-LOCATED HIRS AND MSU MEASURE- 
MENTS PRODUCED BY INGEST* INTO DATASETS FOR DISPLAY AND RETRIEVAL. 
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TOVPRE 
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INPUT (S > 

INGEST OUTPUT < 20 > 
HIRS RTC ( 13 > 

HIRS LCC ( 15 ) 

MSU LCC ( 16 > 

LO-RES T OPOG (17) 
DATA CARD 

IMAGER F ILE ( ?2 ) 
DATA CARD - 
SPECIFY PARAM(S), 
STARTING LINE 


OUTPUT (S ) 

IMAGER F ILE < 22 > 
SOUNDER F I L E < 2 3 ) 
PRINTOUT 


PRINTOUT OF DATA 
COORDINATES 


IK < LINE *ELEMENT> 
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SPECIFY 


WHETHER DATA ARE TO eE ‘LIMB-CORRECTED* CR NOT (SEE BELOW) 


THE PREPROCESSOR PERFORMS THE FOLLOWING FUNC 
IF THE MSU ‘LIMB-CORRECTION* FLAG IS 

ANTENNA PATTERN (TRANSFORM ANTENNA TEH 
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SURFACE REFLECTIVITY (NORMALIZE TO SFC 
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IF THE HIRS ’LIMB-CORRECTION* FLAG IS ON* 
LIMB EFFECTS? 

WATER VAPOR ATTENUATION IN 
IN ADDITION, HIRS CHANNELS 17 
FLUORESCENCE AND REFLECTED 
THE STATE OF THE LIMB-CORP.ECTION FLAG. 
MSU AND HIRS ARE COLOCATED PY INTERPOL AT I 
THE HIRS SCAN PATTERN. 

OUTPUT FILE 22 HAS ALL DATA FOR ONE PARAM 
AND THUS IS OPTIMIZED FOR IMAGING; 
OUTPUT FILE 23 HAS ALL DATA FOP ONE SCAN 
AND THUS IS OPTIMIZED FOR SOUNDING. 

* * .• * • 

* • » t 


THE WINDOW 
AND 18 ARE 
SUNLIGHT. R 


tions: 

MSU DATA ARE CORRECTED FOR 
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HIRS DATA ARE CORRECTED FOR 
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CORRECTED, IN DAYLIGHT* FOR 
ESPECTIVELY, REGARDLESS OF 

NG THE MSU OBSERVATIONS TO 

ETER CONTIGUOUS ON DISK, 

SPOT CONTIGUOUS ON DISK, 


RETRIEVAL 
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FUNCTION: DETERMINE, 

PROFILES OF ATMOSPHERIC 
AS I’ELL AS TOTAL OZONE AND STABILITY 
TIO*’. SEE NOTE ON ’SURFACE DATA* AT 


FROM PREPROCESSED HIRS AND M.SU DATA, VERTICAL 
TEMPERATURE, HUMIDITY, AND GEOPOTENTIAL HEIGHT, 
PARAMETERS, AT HIGH SPATIAL RESQLU- 
THE END OF THIS DOCUMENT. THERE 
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notes: 
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WITH •HRTOPO* . 
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AL FILE* 9 Y LEAST-SQUARES OBJECTIVE ANALYSIS OF HEIGHT FIELDS. 

PROGRAM INPUT (S ) OUTPUT(S) 

;-#*■** ****•«. c * *»** »■<!*». 

WIN'RET RTVL FILE (29) RTVL FILEC29) 

'DATA CARD 

mote: 

TO CONTROL QUANTITY OF PRINTOUT 


b C f % V 

. y * -v 


PEDUNCANCY ELIMINATION 

FUNCTION: ELIMINATE REDUNDANT INFRARED RETRIEVALS* BASED ON VARIABIL 
IT Y IN SELECTED HIRS CHANNELS. 
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rj . , . ■ . .14 

jfc 

Ap: fi ONCXPDRT j 

365 PROGRAM I NPUT (SI GUTPUT(S) ORIGfNAL PAG* Ig 

3/r6 ***•*« >••••*«, Of; POOR QUAU'C/; 

j] |7 RE TRET RT VC FILE(2A> RTVL FILE (29) 

368 DATA CARO 

W 

3ji'c notf: 

3 VI TO SPECIFY OTHER THAN DEFAULT (BUILT-IN) CONTROL PARAMETERS 

372 

Zf^S tms.it- 

Sv’iR (rif M 

375 RETRIEVAL PLOTTERS 

3[,7 FUNCTION: TO PLOT VARIOUS QUANTITIES FROM THE RETRIEVAL FILE. INPUT 

378 CONSISTS OF THE RETRIEVAL FILE(29> AND DATA CAROS TO CONTROL THE LOCA- 
3M T 1 ON AND PARAMETERS PLOTTED ... SEE SOURCE CODE FOR CETAILS. 

3 |0 PROGRAM PRODUCT 

3 el 

3A2 TOVPLFIL 

3 j 13 

3lM TOVPLDIF 

385 TOVPLTEM 

3 Hi, TOVPLDEW 

3 1.-17 TOVPLTBB 

388 

3fp ‘TOVPLFIL* SHOULD EE RUN BEFORE THE NEXT PROGRAM TO BE DESCRIBED 
3jjp (*COMRET*>; THE OTHER PLOTTERS SHOULD BE RUN AFTER FILE-COMPRESSION. 

391 

3 -- f -' r r. 

sljs 

3V9 FILE COMPRESSION 
39J5 

3 j|6 FUNCTION: COMPRESS THE RETRIEVAL FILE BY DELETING SOUNDINGS THAT 

i-7 HAVE BEEN FLAGGED BY *F ILRET* ANO/OR * REDRET* * AND MOVING THE REMAINING 
398 SOUNDINGS TO REPLACE THE ’EMPTY* RECORDS. 


IR-MU SOUNDING DIFFERENCE S, WITH CHARACTERS APPENDED 
TO DENOTE RESULTS OF *FILRET* AND ’REDRET* 

IR-KU DIFFERENCES 

IR-RETRIEVAL TEMPERATURES 

I R-PETR IE VAL DEWPOINTS 

SOUNDING TBB* S FOR SPECIFIED CHANNELS 


PROGRAM INPUT CS ) OUTPUT IS) 

*«.*«*• *«•«•««• 

COMRET RTVL FILE (29 ) RTVL FILE(29) 

**,n 

NOTE ON ’SURFACE DATA* 

PROGRAMS ’TOVRET* AND *ENHRET* INVOKE SUBROUTINE *RTV*, WHICH IN 
TURN CALLS SUBROUTINE * GETSF C * . IN LIKE MANNER, THE ROUTINE *GETSFX* 
IS CALLED FROM WITHIN *FXTIRO*. IT IS OBVIOUS FROM INSPECTION OF THE 
SOURCE CODE THAT THESE ARE DUMMY ROUTINES. THE USER SHOULD PROVIDE AN 
INTERFACE TO ACTUAL SURFACE DATA (1000MB HEIGHT, TEMFERATURE* AND DEW- 
POINT) IF SUCH INFORMATION IS AVAILABLE AT HIS INSTALLATION. 

END OF DOCUMENT 
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TEMPERATURE RETRIEVALS OVER MOUNTAINS USING ONE FOV 
F, Loechner 

Oberpfaf fenhofen DFVLR 
Federal Republic of Germany 

1, Introduction 


Moot of the TOVS software uses boxes of 3x3 fields of view (FOVs) to evaluate 
dear-column radiances. This restricts the horizontal resolution of the 
results to 100x70 km in the subsatellite point. If we can reduce this 
resolution to one FOV, TOVS retrievals would be an important tool in mesoscale 
investigations. Therefore temperature profiles were calculated using 3x3 
boxes and one FOV at sites with different terrain, and compared with 
conventional data. This work was done at the Cooperative Institute for 
Meteorological Satellite Studies (CIMSS) in Madison together with C, M. Hayden 
and A. J. Schreiner. 

2. Technique 

Single FOV retrievals are possible only in cloud free situations. For 5 March 
1982 the AVHRR picture of channel 2 (Figure 1) shows clear sky over parts of 
western Europe and snovcovered Alps. This is the overpass of 1330 GMT, the 
only one, for which DFVLR received AVHRR data. Preliminary comparisons 
between single FOV and 3x3 retrievals over the North Sea and France showed 
excellent agreement. 

The single FOV retrievals were calculated with the program MXTI, the 3x3 
retrievals with the related software MYTI. Both were developed at CIMSS. 

They use the iterative scheme and topographical data of 10 nautical miles 
horizontal resolution. A regression without conventional surface data was 
chosen as a first guess. The main results do not change, if other first guess 
or surface data are applied. 

For the comparisons in mountains, the following groups of adjacent FOV's were 
chosen (Figure 2) : 

a. Lac Leman, an easy terrain without high mountains, with a weather station 
nearby at Payerne and a contiguous FOV in the mountains, 

b. three sites with high mountains at Pontresina, Umbrail pass, and 
Fluchthorn , 

c. four sites around Innsbruck, two of then at the border of the mountains, 
which can be regarded as mixed terrain. Weather stations are at 
Hohenpeissenberg and Garmisch-Partenkirchen. 

3. R esults 

The comparison between the radiosonde soundings of Payerne (Figure 3) and the 
retrievals using single FOV and 3x3 boxes at Lac Leman (Figure A) shows good 
agreement. Details at 550 mb and 300 mb are not reproduced by the 
computations as expected. The deviation of the temperature at 850 mb is 
smaller for the single FOV retrieval. 
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The adjacent FOV at Sion la located in high mountains. In apite of the 
different terrain, the single FOV results for Lac Leman and Sion (Figure 5) 
are very similar. 

There are no radiosonde data for Urabrail pass, Pontresina, and Fluchthorn. 
Assuming that 3x3 retrievals produce correct results, if the type of the 
terrain is the same for all FOV's, only small deviations are found for single 
FOV retrievals (Figure 6) . 

In the comparison between adjacent FOV's with high mountains, Umbrall and 
Pontresina (Figure 7), Pontresina and Fluchthorn (Figure 8), good agreement is 
obtained for single FOV computations. 

Problems arise at the border of the mountains with 3x3 retrievals. Large 
differences are seen between the profiles of Lac Leman and Sion (Figure 9). 

They were not found for the single FOV calculations and ore due to the results 
at Sion. Here the retrievals for Bingle FOV and 3x3 boxes disagree again 
(Figure 10), The comparison with the measurements at Payerne favour the single 
FOV result. 

A similar but less severe situation emerge at Mittenwald for single FOV and 
3x3 retrievals (Figure 11). The observations of Hohenpeissenberg (Figure 12) 
prefer the single TV* profile. Again almost no deviations were found, when 
the single FOV pro/Ves of Mittenwald and the adjacent site Zugspitze are 
compared (Figure 13). 

The reason for the differing results with 3x3 computations seems to be the 
average value of the elevation of the terrain. If the elevations assigned to 
the FOV's used to calculate a 3x3 retrieval differ considerably, problems 
arise. In contrast, retrievals using one FOV produce consistent temperature 
profiles over snow-covered mountains, and at the border of mountains as well. 
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Figure 1. AVHRR picture, Channel ? S March 198 









Fields of view selected for comparisons in mountains 
The diameter of the circles is approximately 20 km. 
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3X3 retrievals at Lac Leman (A) and Sion (B), Figure 10. Single FOV (A) and 3X3 (B) retrieve 

solid lines are temperatures. Sion, solid lines are temperatures. 
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6MT ID TOTRL EQUIP FMfiX CVT LL KXKX PH 
81300 130001 8 0 0 0 -12 4 099399 1 

1151300 130004 A 0 0 0 -13 4 966609 1 

Figure 13. Single FOV retrievals at Mittenwald (A) , 
and Zugspitze (B), solid lines are 
tempera tures. 
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Abstract 


To date, operational satellite temperature retrievals from the TIROS-N/NOAA 
A-C aeries of satellites and a large percentage of those produced for research 
purposes have used statistical techniques to estimate limb effects In satellite 
observed radiances. In this study temperature profiles have been derived using 
the radiative transfer equation in a form which properly takes Into account the 
angle of observation. These temperature profiles have then been compared to 
these derived using the radiative transfer equation with "nadir radiances" 
produced by a statistical limb correction technique similar to those now used 
operationally. This comparison has shown significant differences In the derived 
temperature profiles at large viewing angles, particularly in the case of strong 
meridional temperature gradients. Overall, the results suggest that for the 
calculation of temperature profiles from non-nadir observations, the more proper 
physical solution Is the preferred procedure for deriving temperature fields. 


1. Introduction 

During the "Firat International TOVS Study Conference" at lgls* Austria In 
1983 , stress was put ou the uae of physically based rather than empirical 
algorithms in future temperature retrieval procedures (Menrel and Lynch* 1983). 
In this spirit and also because of the widespread operational use of 
statistically limb corrected data* a study was undertaken to gauge the 
differences between temperature soundings generated from a full physical 
solution of the radiative transfer equation taking into account slant path* as 
opposed to those derived using a statistical correction to the slant path 
radiances followed by a physical solution of a nadir form of the radiative 
transfer equation. 

Two adjacent orbits starting at 1200 GMT and 1340 GMT on 4 March 1982 were 
used for the study. These were chosen as they had been extensively studied at 
the "First International TOVS Study Conference*" overlaid an area with a good 
conventional data base and* most Importantly, exhibited a strong meridional 
temperature gradient. This strong gradient provided a means of examining any 
loss of Information as a result of the limb modelling procedures. 

2. The Retrieval Technique 

The retrieval technique used produces temperature* moisture and ozone 
profiles which are a full solution of the radiative transfer equation. The 
solution methodology is described in detail in Smith et al. (1983). Briefly, 
the physical retrieval algorithm solves for a departure of the true atmospheric 
temperature, moisture* and ozone profile from a "guess" condition, as a function 
of the departure of observed radiances from those calculated for the "guess" 
profile. In this study a climatological first guess was used* tied to a surface 
air temperature and Surface pressure derived from conventional surface data, to 
constrain the lowest levels of the soundings. 


3. The llab Correction Technique* 

3.1 Statletlcel 

The statistical limb correction technique la alnllar to that used 
operationally (Verbowetski, 1981) and vae of the fora used in the "First TO VS 
Export Package." The statistical Hob correction coefficients, vhlch are 
applied to non-nadir radiance observations to give nadir radiances, were 
generated using synthetic radiance data In a manner described by Smith et al. 
(1974). A global representation of 100 temperature and water vapor mixing ratio 
profiles was used to generate synthetic radiances and brightness temperatures at 
nadir and non-nadir angles. In the case of High-resolution Infra-Red Sounder 
(HIRS) observations, each profile was used with a variety of cloud heights and 
amounts to produce the radiances. For the Microwave Sounding Unit (MSU) data, 
although land and water surfaces are treated separately, hybrid samples were 
used for every fifth profile and one In ten soundings over land had randomly 
introduced precipitating clouds. In addition, the effect of sldelobes and 
emisslvity was removed from the nadir MSU radiances. Limb brightness 
temperature correction coefficients were then computed using stepwise multiple 
linear regression, with the differences between the vertical path and slant path 
synthetic data being predlctands and the product of the synthetic non-nadir 
channel brightness temperature and secant of the local zenith angle as well as 
the secant Itself being predictors. 

3.2 Physical 

The radiative transfer equation was used with non-limb corrected data to 
generate temperature profiles by utilizing a computationally efficient 
transmittance model which provided transmit tances appropriate for the 
atmospheric path at the angle of observation. The model algorithm is similar to 
that described by McMlllin and Fleming (1976). For the MSU tran6mittances, the 


effect* of aide lobes on the tranamittances has been computed by e regresaion 
technique rather than using the computationally costly process of convolving the 
antenna gain pattern vlth transmit tances for a number of angles (Smith et al.» 
1983). 

4. The Retrievals 

As previously described, the test orbits were processed using statistically 
limb corrected radiances and non-nadir radiances in conjunction with 
transmit tanccs appropriate to a slant path. The temperatures vere then analysed 
at standard levels using the Barnes analysis scheme (Barnes, 1973). Figure 1 
shows the analyzed 500 mb temperature fields, their differences, and the 
corresponding archived field from ECMWF. It can be seen the thermal gradients 
in the satellite only and ECHWF analyses show good agreement, although a bias of 
1.5 degrees at 500 mb Is evident for these orbits (Smith et al. , 1983). An 
examination of the western edge of the 1340 GMT orbit shows the strong thermal 
gradient region has been best depicted by the full physical scheme. The 
differences between the two satellite analyses show little deviation In the 
center and become quite significant at larger angles. The differences at 700 mb 
and 250 mb are shown in Figure 2. Again, the differences Increase at the edge 
and the large differences appear more frequently as one moves from 250 mb to 
700 mb. The satellite analyses also show a difference over northern Africa 
where we believe the statistical limb correction to be les6 reliable because of 
extreme surface conditions. The distribution of differences was found to be 
similar for the eastern orbit. 

The analyzed temperature fields were compared to the 1200 GMT radiosonde 
data. The differences In the root mean square error and the standard error 


between the llab corrected (LC) end etetietlcelly non- linb corrected (KLC) 
fields (LC-KLC) ere shown in figure 3. In general, both the root seen square 
error end standard error are reduced by using the full physical scheme. At the 
lowest levels the differences are small, partially aa a result of surface 
temperature and pressure constraints on the at those levels, while for the rest 
of the levels a positive lapact is seen as a result of the full physical 
retrieval scheme. 

5. Conclusions 

It has been shown in this study that there appear to be systematic 
differences between soundings generated from statistically limb corrected data 
using a BCheme similar to that now used operationally and those generated using 
transmittance functions appropriate to the angle of view. These differences can 
be significant at larger viewing angles. It appears from the sample examined, 
physical retrievsls which account for the angle of view explicitly are core 
realistic and appear to better maintain temperature gradients in the limb. 
Overall, it would appear that the physical modelling of limb effects in the 
transmittance terms of the radiative transfer equation is the preferred method 
for temperature retrieval. 
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Figure Cep t lone 


Figure i: 


• ) 500 mb temperature anelyele at 1344 GMT on 4 March 1962 

derived uaing observed radiances (NLC) and the full physical 
retrieval technique. 

b) 500 nb temperature analysis at 1344 GMT on 4 March 1982 
derived using statistically limb corrected radiances (LC). 

c) 500 nb temperature analysis from ECMWF for 1200 CMT on 
4 March 1982. 

d) Differences in 500 nb temperature between non-limb corrected 
(NLC) and statistically limb corrected (LC) retrievals. 


Figure 2: a) Differences in 700 nb temperature between non-linb corrected 

(NLC) and statistically limb corrected (LC) retrievals, 
b) Differences in 250 nb temperature between non-limb corrected 
(NLC) and statistically limb corrected (LC) retrievals. 


Figure 3: The differences in root mean square and standard error compared 

to radiosondes between statistically limb corrected (LC) 
retrievals and the non-limb corrected (NLC). 


2/JFL1/12 



FIGURE 1 














Appendix D; The Simultaneous Retrieval Export Package 
W. L. Smith et al. 


excerpted from 

The Technical Proceedings of the Second International 
TOVS Study Conference (editor, P. Menzel) 


THE SIMULTANEOUS RETRIEVAL EXPORT PACKAGE 
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1 . INTRODUCTION 

As part of the "First International TOVS Study Conference", a descrip- 
tion of a physical algorithm for retrieving temperature and moisture sound- 
ings from TOVS radiance data was described (Smith et al., 1983). This 
algorithm was incorporated into the "TOVS Export Package" made available by 
CIMSS to direct readout users of TOVS data. The physical algorithm con- 
sisted of the application of the Smith-Iterative Solution (Smith, 1970) for 
temperature and moisture profiles. In that approach, (1) an initial profile 
for temperature and water vapor is specified from climatology, statistical 
regression, or from an KWP model, (2) radiances are calculated from the 
initial profiles, (3) the temperature profile is adjusted in on iterative 
manner until there Is agreement between the observed and calculated radi- 
ances in the cloud insensitive microwave 0* channels, (4) the infrared 
window channels are used to define either the surface skin temperature or 
the temperature of cloud within the instrument's field of view and the cloud 
level is defined using the microwave specified temperature profile, (5) the 
guess moisture profile is adjusted to reflect the existence of cloud by 
assuming 100? relative humidity at the cloud level and then further adjusted 
in an iterative manner in order to achieve convergence between observed and 
calculated radiance for the water vapor channels, and (6) the temperature 
profile is then further adjusted in an iterative manner in order to achieve 
convergence between the radiances observed and calculated in the infrared 
CO^ channels. 

Since the first Xgls conference, a new retrieval algorithm has been 
developed which permits the simultaneous retrieval of surface-skin (or 
cloud) temperature, and the temperature and moisture profiles. The 
advantage of the "simultaneous solution" is two-fold: (1) the radiances 

observed in all channels are used to solve for all parameters simultaneous- 
ly, thus alleviating the problem of the interdependence of the radiance 
observation upon the three parameters, and (2) since a direct analytical 
solution is employed, the process is computationally efficient. 

In this paper the new algorithm now incorporated in the "TOVS export 
package" is described. As with the prior "iterative" algorithm, the 
physical nature of the solution permits the influence of surface variables 
(i.e., terrain elevation, emissivity, and temperature) and cloudiness to be 
accounted for in the profile determinations. The cloud handling algorithm 


is modified to enable the Infrared data to be utilized in partly cloudy as 
well as cloud overcast sky conditions. The low sensitivity of the simultan- 
eous solution to the initial guess profile and the improved performance of 
the simultaneous retrieval method compared to the previously established 
iterative technique (Smith et. al., 1983; Sussklnd, 1984) are demonstrated 
from TOVS orbits over Europe obtained during the ALPine Experiment (ALPEX) 
as selected by the International Radiation Commission's TOVS Working Group 
for the lntercomparlson of retrieval methods. The physics for treating 
surface emlssivity, terrain elevation, and reflected sunlight are not 
described here since these aspects have been provided in the previous report 
(Smith et al., 1963). 


2. DIRECT PHYSICAL SOLUTION 


An important advance in the profile retrieval methodology is the simul- 
taneous temperature and water vapor solution (Smith and Woolf, 1984). In 
order to achieve the simultaneous solution, the integral form of the 
radiative transfer equation is integrated by parts and treated in the 
perturbation form: 


3T 8t (3B/3T) , S ^ il ( <*E/ <*T) ^ 

6T " / 61 \ »U ( 3B/9T*) dp ‘ / ^Sp (W9T*) dp (3B/9T*) (1) 


A 

where T is brightness temperature, U is precipitable water vapor, B is 
Planck radiance, T is temperature, T is surface-skin temperature, t is 
transmittance, and p is pressure. T$e perturbation, 5, is with respect to 
an a-priori estimated or mean condition. The pressure dependence of all 
integrand variables is to be understood. In order to solve (1) for 4U, 4T, 
and 6 t from a set of spectrally independent radiance observations, the 
perturlation profiles are represented in terms of arbitrary pressure func- 
tions, $(p) : 
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6T(p) - - E « ♦ (p) (2b) 
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where q(p) is the water vapor mixing ratio and g is gravity. The zero sub- 
script indicates the a priori condition. Equation 2a implies from the gas 
lav: and hydrostatic equation that 
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Substituting representations (2b) and (2c) into (1} and letting ap=6T 
yields for each spectral radiance observation, fiT , for a set of K spectral 
channels: 
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The 4> . quantities are calculated from the a- priori estimated or mean 

prof il^^condltions. Written In matrix form 

t* - $a (5) 

* 

where t is a row vector of K radiance observations, a is a row vector of 
M+l coefficients, and $ is a matrix having dimensions (K>4i+1). Assuming 
that K2H+1, then the least squares solution of (5) is employed to give 


t-1t*'vT -IT* 

a m ($*$) t t * + yl) * t (6) 

T -1 

where ( ) indicates matrix transposition and ( ) indicates matrix 

inverse. The yl term, where y (nominally 0.1) is a scalar and 1 is the 

identity matrix, is incorporated to stabilize the matrix inverse. Once the 

* 6 are determined, 6T , fiq, and 6T are specified from (2) and added to the 

a-priori estimates to yield the final solutions for surfBce-skin temperature 

and the water vapor mixing ratio and temperature profiles. 

The choice of pressure basis functions $(p) is arbitrary. For example, 
empirical orthogonal functions (i.e.» eigenvectors of the water vapor and 
temperature profile covariance matrices) can be used in order to include 
statistical information in the solution. However, here for its application 
to the TIROS Operational Vertical Sounder (TOVS) data ana the VISSR 
Atmospheric Sounder (VAS) observations, physical functions, the profile 
weighting functions (dt/dlnp) of the radiative transfer equation, are used 
as the basis functions. 

Ancillary information such as surface observations can be easily 
incorporated into the profile solutions. For example for surface observa- 
tions it follows from (2) that 
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have the same form as (3) and therefore can be added to the aet to yield K+2 
equations to solve for M+l unknowns (a). 


The main advantage of the "simultaneous" retrieval method is that it 
enables the temperature and water vapor profiles and the surface skin 
temperature to be determined simultaneously using all the radiance informa- 
tion available. This solution alleviates the problem associated with water 
vapor radiance dependence upon temperature and the dependence of several of 
the carbon dioxide channel radiance observations used for temperature 
profiling, on the water vapor concentration. The dependence of the radiance 
observations on surface emissions is accounted for in the simultaneous 
solution by the inclusion of surface temperature as an unknown. Also, since 
only a single matrix inversion is required for the specification of all 
parameters, the solution is computationally efficient. Finally, ancillary 
observations of temperature and/or moisture from surface sensors or air- 
craft, for example, can be readily Incorporated in the solution. 


0 

n 
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Logistics 

The processing of the TOVS data begins with the specification of an 
initial profile and surface observations, if available. The initial profile 
can be specified either from climatology or through the use of regression 
coefficients, based either on synthetic radiances or on real radiances 
matched to radiosonde observations. The climatological and regression 
generated initial profile algorithms are Internal to the Export Package. A 
third option exists whereby an analysis or forecast generated grid field of 
temperature and water vapor profiles can be tased to create Initial profiles. 
In this case the grid fields are external to the Export Package as are the 
surface data grid fields. 


H Once a guess profile is generated, MSU and HIRS radiances are obtained 

[ ] for a particular 3x3 array of FIRS FOV's, as described in the. earlier 

report (Smith et al., 1983), the MSU data having been interpolated to each 
r~> HIRS FOV. As described in the next section, a search through the nine FOV’s 

1 1 is performed to arrive at a set of "clearest" radiances for the sounding 

1 retrieval. 


The stratospheric HIRS brightness temperatures (channels 1-3), the MSU 
brightness temperatures (MSU 1-4), and the middle and upper tropospheric 
HIRS water vapor channels (HIRS 11 and 12) are used to derive a first 
estimate of the temperature and moisture profile for the sounding location. 
Because of the channels used, this estimate should be relatively free of 
error due to cloud attenuation. The weighting functions for HIRS-1 and MSU 
2-4 are utilized as the basis functions for temperature and the HIRS 7 and 


12 weighting functions are used as the basis functlonas for water vapor in 
the Initial retrieval. Once the first estimate of the temperature and water 
vapor profile Is achieved, the height and amount of any cloud affecting the 
Infrared observations Is determined by the method described In the next 
section. After this Is accomplished all channels (except H1RS channels 
13-19 in the cloudy case) are used to calculate the final surface tempera- 
ture and the temperature and water vapor profiles. For the achievement of 
the final profile estimates the weighting functions for HIRS 1, 3, 7 and HSU 
2-4 are used ss the temperature profile basis functions and those for HIRS 
7, 11, and 12 are used as the water vapor profile basis functions. Since as 
many ss nineteen different spectral radiance observations are used for the 
surface temperature and profile retrievals, the system of equations to be 
Inverted is heavily overdetermlned, thereby stabilizing the solution. (It 
should be noted that ozone and geopotential height are determined in the 
came manner as described for the iterative solution Export Package.) 


Handling the Influence of Clouds 


In the TOVS data processing, soundings are derived from a 3 x 3 matrix 
of HIRS fields of view (FOV), the Microwave Sounding Unit (MSU) data being 
spatially Interpolated to the location of the HIRS spots. From the array, 
the observations for the "clearest" FOV's, defined as those whose 1 1 pm 
radiance values are within 2°K of the maximum, are averaged for the sounding 
determination. The magnitude of the visible channel reflectance and the 
3.7 vm, 4.0un, and 1 1 pan window channel brightness temperatures are used in 
conjunction with surface temperature observations, if available, to specify 
whether the "clearest" radiances are contaminated by clouds. If the radi- 
ances are determined to be affected by clouds then their pressure height and 
fractional coverage are specified on the basis of certain CO^ channel 
infrared radiances and the microwave radiance specified temperature profile. 
The cloud height is calculated using the GO- slicing technique (Smith and 
Platt, 1977, Menzel et al., 1983) from the relation 
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where I (2) is the clear-column .radiance calculated from the microwave 
specified temperature profile, T, and v and v refer to HIRS channels 5 and 
7 respectively. The cloud pressure P is obtained by trial and error. The 
cloud fraction, N, is then obtained from the relation 


(Vj) is Che radiance calculated for an opoque cloud overcast 

a 

Given P and N, the guess mixing ratio profile is adjusted by assuming 
that the mixing ratio at the cloud level Is given by 

W ( P C ) “ NW Sflt fT(p c >) ♦ (1 - N)W(p c ) (9) 

where W is the saturation mixing ratio corresponding to the microwave 
specif iea c temperature at the cloud pressure P and U 1 b the original guess 
value of mixing ratio. Belov the cloud a nev c guess mixing profile Is 
achieved by Interpolation using the original surface mixing ratio value. 

Once the cloud parameters and the guess mixing ratio is established 
then the brightness temperature discrepancies from the guess conditions, 
dT*, can be calculated, including the effects of cloud within the field of 
view of the H1RS instrument. Also, for the infrared channels the effects of 
cloud and the transmittance function, t, employed in equation (4) can be 
properly included by assuming that 

t - (1 - N) t or * 8 

for pressures greater than the cloud pressure. This latter transmittance 
function modification is important to reduce the Influence of the infrared 
observations on the solution below the cloud in pxoportion to the cloud 
obscuration. Once this step is employed, the simultaneous solution proceeds 
in exactly the same manner as in the clear sky condition. 
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